Fast nucleon emission as a probe of the isospin momentum dependence 
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In this article we investigate the structure of the non-local part of the symmetry term, that 
leads to a splitting of the effective masses of protons and neutrons in asymmetric matter. Based 
on microscopic transport simulations we suggest some rather sensitive observables in collisions of 
neutron-rich (unstable) ions at intermediate (RIA) energies. In particular we focus the attention on 
pre-equilibrium nucleon emissions. We discuss interesting correlations between the N/Z content of 
the fast emitted particles and their rapidity or transverse momentum, that show a nice dependence 
on the prescription used for the effective mass splitting. 
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I. INTRODUCTION 

Effective interactions have been extensively used in 
transport codes to gain knowledge about nuclear mat- 
ter properties in conditions far from equilibrium. They 
predict very different behaviours for some physical pro- 
perties, such as the momentum dependence of the mean 
field in the iso-scalar and iso- vector channel (see Q, Q for 
recent reviews). 

In general, the momentum dependence of the potential 
is strictly related to the concept of in-medium reduction 
of the nucleon mass. In particular, the non-local part of 
the symmetry potential has become an interesting sub- 
ject of investigation, because it leads to a splitting of the 
masses of different nucleonic species in asymmetric mat- 
ter, and the sign of this splitting already is an open and 
controversial problem (0,0,11). 

Considering separately the local and non-local contri- 
butions to the potential part of the symmetry energy, 
a sharp change is observed, for instance, going from 
the earlier Skyrme forces to the Lyon parametrizations, 
with almost an inversion of the signs of the two con- 
tributions Q{. The important repulsive non-local part 
of the Lyon forces leads to a completely different be- 
havior of the neutron matter EOS, of great relevance 
for the neutron star properties. Actually this substan- 
tially modified parametrization was mainly motivated by 
a very unpleasant feature in the spin channel of the earlier 
Skyrme forces, the collapse of polarized neutron matter 
0,0, 00- In correspondence the predictions on isospin 
effects on the momentum dependence of the symmetry 
term are quite different. A very important consequence 
for the reaction dynamics is the expected inversion of the 
sign of the n/p effective mass splitting, i.e. in the Lyon 
forces neutron effective masses are below the proton ones 
for n-rich matter. 

We note that the same is predicted from microscopic 
relativistic Dirac-Brueckner-Hartree-Fock (DBHF) cal- 



culations 0, 0, and in general from the introduc- 
tion of scalar isovector virtual mesons in Relativistic 
Mean Field (RMF) approaches [H At variance, 
non-relativistic Brueckner-Hartree-Fock (BHF) calcula- 
tions are leading to opposite conclusions 0, 0, Il7j . 
However, the comparison between relativistic effective 
(Dirac) masses and non-relativistic effective masses re- 
quires some attention, see the Ch.6 of 2] and refs. 
therein. In fact, though the neutron Dirac-mass is be- 
low the proton one, in correspondence the non-relativistic 
Schroedinger-m&ss splitting can have any sign. Indeed 
this is just the case of the very recent DBHF analysis 
of the Tubingen group 0, 0|. Such puzzling effect is 
related to the intrinsic momentum dependence of the nu- 
cleon self-energies, it is then beyond the RMF approach 
and it represents a very sensitive test of the treatment of 
Dirac-Brueckner correlations. All that clearly shows that 
sensitive experimental observables are largely needed. 

The sign of the splitting will directly affect the energy 
dependence of the so-called Lane Potential, i.e. the dif- 
ference between (n,p) optical potentials on charge asym- 
metric targets, normalized by the target asymmetry [ljj. 
A decreasing behaviour is observed in the case m* > m*, 
while a positive slope is obtained in the opposite case. 
An important physical consequence of the negative slopes 
is that the isospin effects on the optical potentials tend 
to disappear at energies just above 100 MeV (or even 
change the sign for "old" Skyrme-like forces). Also, we 
expect a crossing of the two prescriptions at low ener- 
gies, i.e. low momentum nucleons will see exactly the 
same Lane potentials, as shown in detail in ref. |5(. 

Unfortunately results derived from neutron/proton 
optical potentials at low energies are not conclusive 
[la. I20I |2"H , since the effects due to the mass splitting 
appear of the same order of the uncertainty on the de- 
termination of the local contribution to the symmetry 
energy. Moreover at low energies, due to the crossing 
discussed before, it is difficult to appreciate differences 
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between positive or negative slopes. 

A positive slope is obtained, for instance, within the 
phenomenological Dirac Optical Potential (Madland — 
potential), with different implicit momemtum depen- 
dences in the self-energies |22J, |23j. This potential has 
been constructed fitting simultaneously proton and neu- 
tron (mostly total cross sections) data for collisions with 
a wide range of nuclei at energies up to 100 MeV. Re- 
cently such Dirac optical potential has been proven to 
reproduce very well the new neutron scattering data on 
20S Pb at 96 MeV measured at the Svendberg Labo- 
ratory in Uppsala. 

More data are needed at higher energies 
(around/above 100 MeV), to improve the system- 
atics, in order to clearly disentangle between the two 
trends of the Lane potential and of the effective mass 
splitting. 

We note however that these properties of the interac- 
tion will also affect the dynamical evolution of heavy ion 
collisions. So we can get indipendent information about 
it just by looking at some suitable reaction observables. 
We can expect important effects on transport proper- 
ties ( fast particle emission, collective flows) of the dense 
and asymmetric NM that will be reached in Radioactive 
Beam collisions at intermediate energies. 

Here we will focus on the study of pre-equilibrium 
emission in neutron-rich central collisions, at 50 and 100 
MeV/A. We will see that the energy dependence of the 
N/Z content of fast emitted particles is particularly sen- 
sitive to the sign of the effective mass splitting. Hence it 
would be possible to answer many fundamental questions 
in isospin physics by looking at appropriate observables 
in intermediate energy HIC. 



(Z), and the total mass A — N + Z; the quantity go is 
the normal density of nuclear matter. The momentum 
dependence is contained in the X TT i terms, which indicate 
integrals of the form: 



2"tt' = / dp dp' f T {r,p)f T >(r,p')g(j),j)') 

with g(p,p') = (p ' — p') 2 ■ This choice of the function 
g{p,p') corresponds to a Skyrme-like behaviour and it is 
suitable for BNV simulations. We use a soft equation 
of state for symmetric nuclear matter (compressibility 
modulus Knm{Qo) = 215 MeV). In this frame we can 
easily adjust the parameters in order to have the same 
density dependence of the symmetry energy but with two 
opposite n/p effective mass splittings, as predicted by 
the early Skyrme forces Q and the later Skyrme-Lyon 
parametrizations Q . So we can separately study the cor- 
responding dynamical effects, 

A good way to visualize the physical meaning of the 
mass splitting is to look at the kinetic energy dependence 
of the Lane potential Uhane — (U n — U P )/2I, 0. Its 
slope depends on the value and sign of the mass splitting 
(see |5| for details). In Fig. Q]we plot the Lane potential 
with the parameters corresponding to the two choices of 
the sign of the n/p mass splitting (shown in the insert for 
the / = 0.2 asymmetry). In the two cases, the absolute 
value of the splitting is exactly the same (the difference 
between neutron and proton masses is ~ 10% at normal 
density), only the sign is opposite. The upper curve well 
reproduces the Skyrme-Lyon (in particular SLyA, Sly7) 
results, the lower (dashed) the SIII, SKM* ones. 



II. DETAILS OF THE INTERACTION 

In order to directly test the influence of the 
Schroedinger-effective mass splitting we will present re- 
sults from reaction dynamics at intermediate energies 
analysed in a non-relativistic transport approach, of 
BNV type, see [lE^. The I so - MD effective interac- 
tion is derived via an asymmetric extension of the GBD 
force IMS • 

The energy density can be parametrized as follows (see 
also 0,13): 

e = e km + e(A', A") + e(B', B") + e(C , C") (1) 
where Skin is the usual kinetic energy density and 

e(A',A") = (A' + A"I 2 ) — 

eo 

s(B',B") = (B' + B"I 2 )(j^j g 
e(C, C") = C'(I NN + Ipp) + C"T NP (2) 

The variable / = (N — Z) / (A) defines the isospin content 
of the system, given the number of neutrons (N), protons 
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FIG. 1: Energy dependence of the Lane potential for the used 
parametrizations. Small panel: the effective mass splitting, 
related to the slope of the Lane potential (for a L = 0.2 asym- 
metry). 

From the figure we can immediately derive the ex- 
pectation of very different symmetry effects for nucle- 
ons around lOOMeF kinetic energy: enhancement (with 
larger neutron repulsion) in the m* < m* case vs. a 
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disappearing (and even larger proton repulsion) in the 
to* > to* choice. Besides, we can see a crossing of the 
two prescriptions at low energy: this means that low mo- 
mentum nucleons feel almost the same Lane potential. 
Thus it is important to choose the appropriate dynamical 
observables in order to clearly sec the effect of mass split- 
ting, e.g. to look at observables where neutron/proton 
mean fields at high momentum are playing an important 
role. 



III. ISOTOPIC CONTENT OF THE FAST 
NUCLEON EMISSION 

We have performed realistic "ab initio" simulations of 
collisions at intermediate energies of n-rich systems, in 
particular 132 Sn + 124 Sn at 50 and 100 AMeV, at im- 
pact parameter b = 2 fm. We have employed either 
an asy — stiff or an asy — soft density dependence of 
the symmetry energy, corresponding to different slopes 
around normal density, with an increasing symmetry re- 
pulsion in the asy — stiff case, |28j . One of the aims 
of our work is also to select reaction observables more 
sensitive to the momentum (non-local part) than to the 
density (local part) dependence of the symmetry term. 

Performing a local low density selection of the test par- 
ticles (p < po /8) we can follow the time evolution of nu- 
cleon emissions (gas phase) and the corresponding asym- 
metry. It is shown in Fig. [2 (where I(t) = (N - Z)/A) 
for the 50 AMeV reaction. The solid line gives the initial 
asymmetry of the total system. 
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FIG. 2: Isospin gas content as a function of time in a central 
collision 132 Sn + 124 Sn at 50 AMeV for two opposite choices 
of mass splitting. Left panel: asy-soft EOS; right panel: asy- 
stiffEOS. 

A general feature of nucleon emission is the isospin 
distillation (see |2^| and refs. therein), which leads to a 
gas phase with an asymmetry greater than the initial one. 



We find an interesting behaviour for the isospin content of 
the nucleonic gas: at early times (up to about 60 fm/c), 
during the (high density) compression phase, the asy- 
stiff choice leads to higher isospin asymmetry than the 
asy-soft one; later, we can see an inversion of this trend, 
due to large contributions to the gas phase coming from 
low density regions [28j; finally , th e gas phase is more 
asymmetric in the asy-soft case ,29j. 

This behaviour is mainly due to the local part of the 
symmetry energy, but for both choices of the density stiff- 
ness of the symmetry term we clearly see also the effects 
of the mass splitting, resulting in a reduced fast neutron 
emission when to* > m*. It is more pronounced during 
the early stages of the reaction (up to 60 fm/c), when the 
most energetic particles are emitted [30j | and the momen- 
tum dependent part of the mean field is more effective. 
The two effects we have just discussed, the former related 
to the density dependence, the latter to the momentum 
dependence of the symmetry potential, can be singled 
out in the transverse momentum distribution of the N/Z 
content of the gas phase. 

In order to better isolate the mass splitting effect and 
to select the corresponding observables, in Figs. |21 El 
we report the N/Z of the "gas" at two different times, 
t = 60 fm/c (end of the pre-equilibrium emission), and 
t = 100 fm/c (roughly freeze-out time). We have followed 
the transverse momentum dependence, for a fixed central 
rapidity y^ (normalized to projectile rapidity). 
In Fig. |3| we show the results without the mass splitting 
effect, i.e. only taking into account the different repulsion 
of the symmetry term (the initial average asymmetry is 
N/Z = 1.56). 
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FIG. 3: Transverse momentum dependence of the neu- 
tron/proton ratio in the rapidity range \y \ < 0.3 for a cen- 
tral reaction 132 Sn+ 124 Sn at 50 AMeV. Comparison between 
two choices of symmetry energy stiffness at the two different 
times defined in the text. 

At early times (left panel of Fig. when the emission 
from high density regions is dominant, we see a differ- 
ence due to the larger neutron repulsion of the asy-stiff 
choice. The effect is reduced at higher p/s due to the 
overall repulsion of the isoscalar momentum dependence. 
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FIG. 4: Same as Fig. /or too opposite choices of mass 
splitting. 



FIG. 5: Same as Fig. ^| but for an incident energy of 
100 AMeV (note the reduced "pre- equilibrium" , left panel, 
and "freeze- out" , right panel, times). 



Finally at freeze out (right panel) the difference is less 
pronounced, since Isospin Distillation is coming into 
play during the expansion phase; it is more efficient in 
the asy-soft choice, and refs. therein. 

When we introduce the mass splitting, the difference 
in the isotopic content of the gas, for the larger trans- 
verse momenta, is very evident at all times, in particular 
at the freeze-out of experimental interest, see Fig. 0] 
As expected the m* > to* sharply reduces the neutron 
emission at high pt's. On the other hand, low momen- 
tum nucleons show a behaviour which is consistent with 
Fig. |21 mainly ruled by the local part of the symmetry 
potential. 

We have repeated the analysis at higher energy, 
lOOAAIeV, and the effect appears nicely enhanced, see 
Fig. [S] The observed differences for nucleons at high p t 's 
already appear at early times: since these momentum 
regions are populated mostly during the early stages of 
the reaction (cfr. [13), we finally get a permanent signal 
of a dynamical feature of the interaction. We note a de- 
creasing behaviour of the N/Z content of pre-equilibrium 
emission versus the transverse momentum, due to a larger 
neutron attraction, in the to* > to* case, while in the op- 
posite case we get a flat, or even increasing, trend. 

The N/Z of emitted particles has been studied also as a 
function of rapidity (see FigHJJ). It is possible to see that, 
at 50 /to / c (end of the pre-equilibrium emission at this 
beam energy), differences between the two prescriptions 
are present at all rapidities, though more pronounced for 
faster particles, as expected. We note the opposite be- 
haviour of the two mass splitting sign. Going to 80 fm/c 
(new freeze-out time), the differences at small rapidity al- 
most disappear, while the different effects already present 
at 50 fm/c for large rapidities are kept and are even en- 
hanced. The behaviour observed as a function of rapidity 
is quite similar to the trend already seen as a function 
of the transverse momentum. The different mass split- 
ting sign leads to different results, expecially for particles 
with larger rapidity or, in a given rapidity bin, with larger 
transverse energy. Using the to* < to* prescription, the 
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FIG. 6: The asymmetry of fast emitted particles as a function 
of rapid; 
b=2 fm. 



of rapidity, for the reaction 132 Sn + 124 Sn at 100 MeV/A, 



N/Z content of fast emitted particles is enhanced, while 
it is reduced in the opposite case. 

This is fully consistent with the expectations deduced 
from the Fig. ^at the end of Sect. II. 

Similar results have been obtained in ref. |3lj | at 
AQQAMeV in a different I so — MD model, restricted 
to the to* > to* choice. 



IV. CONCLUSIONS 

Our aim has been to find an observable effect of the 
nucleon effective mass splitting, which is determined by 
the momentum dependence of the symmetry potential. 
We have performed transport simulations for reactions 
of interest for the new radioactive beam facilities at in- 
termediate energies. 

From our results it appears that the isospin content of 
pre-equilibrium emitted nucleons at high transverse mo- 
mentum, or large rapidity, is rather sensitive to this prop- 
erty of the mean field. Thus it can be a good observable 
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in order to learn about fundamental properties of the nu- 
clear interaction, such as its momentum dependence in 
the isovector channel. In particular, it is found that the 
m* > m* splitting leads to a decreasing behaviour of the 
N/Z of emitted particles versus rapidity or transverse en- 
ergy, while the opposite splitting is associated with a flat 
(or even increasing) trend. 



We finally like to remark that this is an observable of ex- 
perimental interest that can be studied, even in absence 
of information about neutron emissions, by looking at the 
correlation between kinematical properties and isotopic 
content of light clusters emitted during the early stage of 
the collision. 
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